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Reconstructed jets in heavy ion collisions are a crucial tool for understanding the quark-gluon 
plasma. The separation of jets from the underlying event is necessary particularly in central heavy 
ion reactions in order to quantify medium modifications of the parton shower and the response of the 
surrounding medium itself. There have been many methods proposed and implemented for studying 
the underlying event substructure in proton-proton and heavy ion collisions. In this paper, we detail 
a method for understanding underlying event contributions in Au+Au collisions at y^ijvjv = 200 
GeV utilizing the HIJING event generator [l|. This method, extended from previous work by the 
ATLAS collaboration [3| , provides a well-defined association of "truth jets" from the fragmentation 
of hard partons with "reconstructed jets" using the anti-fcr algorithm. Results presented here are 
based on an analysis of 750M minimum bias HIJING events. We find that there is a substantial 
range of jet energies and radius parameters where jets are well separated from the background 
fluctuations (often termed "fake jets") that make jet measurements at RHIC a compelling physics 
program. 



I. INTRODUCTION 

Understanding the detailed interaction and coupling 
between hard scattered partons and the quark-gluon 
plasma through which they propagate is essential to our 
fundamental knowledge of QCD and in determining prop- 
erties of the quark-gluon plasma. The measurement of 
fully reconstructed jets in heavy ion collisions at the 
LHC [3, 4] highlight the substantial additional informa- 
tion contained therein and its complementary nature to 
single hadron [5—7] , di-hadron correlations [§-Hlj. The 
measurement of direct photon-jet correlations is another 
critical handle to be utilized [12| . Extending fully calori- 
metric jet measurements to lower center-of-mass energies 
at the Relativistic Heavy Ion Collider provides measure- 
ments for kinematics difficult to access at the LHC and 
the QGP at different temperature and coupling regime. 

With the first Pb+Pb at s /sWn=2. 76 TeV collisions at 
the LHC new insights into jet physics in heavy ion col- 
lisions were gained. The ATLAS collaboration reported 
an increase in the number of energy asymmetric di-jets 
in central Pb+Pb collisions compared to proton-proton 
and peripheral Pb+Pb collisions [3j . They also reported 
the suppression of jets with 100< pr <200 GeV/c by a 
factor of approximately when comparing central to pe- 
ripheral Pb+Pb collisions [ljj. The CMS collaboration 
measured jet-hadron correlations in a similar jet pt range 
and found that the energy lost by high pt fragments was 
approximately balanced by very low px tracks far from 
the jet axis [4J. However the data from both RHIC and 
the initial LHC results are not enough to constrain the 
physics of jet quenching. Most theoretical descriptions 
have relied on weakly coupled techniques [14j. Features 



of strong coupling, as observed in descriptions of the bulk 
matter, might contribute to jet quenching as well. More 
data on jet observables (including dijet, 7-jet and heavy 
flavor tagged jets) at RHIC and the LHC will be neces- 
sary to understand the physics of jet quenching over the 
full range of medium properties and jet kinematics and 
probe for sensitivity of the quenching to outgoing parton 
virtuality. 

The multiplicity of charged particles dN c h/drj is ap- 
proximately 2.15 times higher for Pb+Pb central colli- 
sions at y/SNN = 2.76 TeV compared with Au+Au cen- 
tral collisions at ^/snn =200 GeV [la ]. Thus the soft par- 
ticle background is substantially higher for LHC events. 
However, the jet cross section is substantially higher as 
well, and measurements for jets with energies greater 
than 100 GeV appear well separated from the background 
(though detailed publications of these studies are not yet 
available). Various methods have been explored at the 
LHC and RHIC for understanding the underlying event 
contributions, and what are often referred to as "fake 

jets" [lEMa. 

At y/sNN=200 GeV the projected jet rates into \rj\ <1 
based on NLO pQCD cross sections I20il an d expected 
RHIC luminosities have been computed [21j. In a typi- 
cal year of RHIC running 50B Au+Au events could be 
sampled. In the top 20% centrality that would lead to 
approximately 10 7 jets above 20 GeV, 10 6 jets above 30 
GeV, 10 5 jets above 40 GeV and 10 4 jets above 50 GeV. 
Over 60% of the time there is full containment of the 
opposing dijet for 20 GeV jets with that percentage in- 
creasing with increasing jet energy. 

In this paper, we present a study of jet reconstruc- 
tion and separation from the underlying event using 



the HIJING [l] event generator for Au+Au events at 
y/sNN=200 GeV. This follows an iterative underlying 
event subtraction procedure extended from one devel- 
oped by the ATLAS Collaboration (2J. While the exact 
definition of a correctly reconstructed jet versus a "fake 
jet" is arbitrary, this methodology allows us to make a 
well-defined and documented comparison to cross-check 
with other methods. 



II. JET - UNDERLYING EVENT SEPARATION 
METHODOLOGY 

For these studies we utilize the HIJING (version 1.383) 
event generator run with standard settings and quench- 
ing turned off for Au+Au collisions at y / JJv r /v=200 GeV. 
HIJING is a QCD based Monte Carlo for the study of jet 
production in high energy nucleus-nucleus collisions. For 
these initial studies, we explore what a "perfect" detec- 
tor is capable of measuring. We assume a segmentation 
in At? x A<j> = 0.1x0.1 and that all particle energies are 
recorded perfectly (with the exception of neutrinos and 
muons). We assume a nominal coverage of \t]\ < 1.0 and 
full azimuthal coverage. For the entire study we utilize 
the anti-fey jet reconstruction algorithm 22] (part of the 
FastJet package (23J) with radius parameters R=0.2, 0.3, 
and 0.4. 

A schematic diagram of the underlying event subtrac- 
tion steps is shown in Fig. [T] The first step is to run 
the anti-fey algorithm over the full set of energy values 
(unsubtracted) with R=0.2 and record the jet axis coor- 
dinates in r\ and cj). This initial suite of jet candidates is 
used to exclude regions around these jets from the ini- 
tial underlying event average energy. Exclusion regions 
are defined by R=0.2 jets in which the maximum tower 
in the jet has an energy of more than three times the 
average tower energy in the jet. We then exclude all en- 
ergies for 0.1 x 0.1 cells whose center coordinate is within 
AR = ^/ Ar\ 2 + Ac/) 2 < 0.4 of any of the above initial 
jet candidates. The remaining energy values are used 
to determine the average cell energy (i.e. in the non- 
jet regions) in A77 = 0.1 strips. The modulation in the 
background due to flow must first be removed, so the 
(cos(2(/))) (i.e. i>2 parameter) is also determined for the 
energy distribution and removed from each cell before 
determining the average. Only the 77 strips which have 
complete <f> coverage after the determination of the ex- 
clusion regions are used in the vi determination. 

The HIJING generator has no bulk collective flow and 
thus has only a modest (cos(2ip)) from decays, di-jet cor- 
relations and global momentum conservation. As the flow 
modulation of the underlying event is an important com- 
ponent of any subtraction procedure on real data, we 
have added a flow modulation to the individual HIJING 
particles prior to segmenting the energies into cells. The 
flow parameterization [24j is based on fits to the available 
data. Higher flow moments have an increasing relative 
importance for more central events [25|, and can be in- 



corporated in future studies. 

This underlying event average energy is a zeroth order 
estimate since the initial jet determination does not have 
an underlying event subtraction. We now subtract the 
V2 modulated underlying event energy cell-by-cell from 
the cells contained by the initial set of R=0.2 jets to get 
a better estimate of the jet Et- At this second iterative 
step, new exclusion regions arc defined by towers with 
AR < 0.4 around background subtracted jets with Et > 
20 GeV. The underlying event and V2 are re-determined 
removing towers within Ar] <0.4 of the jets as described 
above and the background re-subtracted from the orig- 
inal unsubtracted towers. Finally the anti-fey jet algo- 
rithm is run on the background subtracted towers with 
a range of R values (0.2, 0.3, 0.4). 

When the jet reconstruction is run over background 
subtracted towers many of the towers have negative Et- 
We modify these towers to have a small positive energy 
before passing them to the jet reconstruction algorithm. 
After the towers are grouped into jets we recalculate the 
jet Et including the negative energy. 



III. HIJING TRUTH INFORMATION 

In order to identify "true jets" from the HIJING event 
generator, we have augmented the code so that every 
time the fragmentation routine (HIJFRG) is called, we 
record the set of final state hadrons that result from that 
fragmentation. We then run the anti-fey algorithm on 
those final state hadrons (using their exact momentum 
vectors). The jet reconstruction is run once for each anti- 
fey jet R parameter under consideration and the resulting 
"true jet" information is added to the output. 

Before presenting the results, it is important to define 
our terms. Even in a model such a HIJING where all 
truth information is known, there is an arbitrariness in 
the definitions of "fake jets" and "true jets" as examples. 
For example, consider a HIJING fragmentation call that 
results in hadrons reconstructed via anti-fey with R = 
0.4 into a jet with an energy of 20 GeV. If after running 
jet reconstruction over the full HIJING event one recon- 
structs a jet using anti-fey with R = 0.4 that has a jet axis 
within AR = ^J Arf + A(f> 2 < 0.25 and energy 18 GeV, 
common sense might dictate that this was a "true jet" 
and the 2 GeV difference is a result of the fluctuations in 
the underlying event. However, imagine that the HIJING 
fragmentation reconstructed to an energy of 4 GeV and 
the full HIJING event resulted in a jet along the same 
axis but with energy of 40 GeV. Common sense might 
dictate that this was a "fake jet" (i.e. a very small jet 
that combined with substantial background fluctuations 
that results in a very large energy reconstruction). 

Here, we define a "fake jet" as one where the associated 
HIJING fragmentation jet is less than 5 GeV (or does 
not exist at all). The jet is a good "true jet" if there 
is an associated HIJING fragmentation jet within AR = 
\J Ajf + A(j> 2 < 0.25 and greater than 5 GeV. We then 
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FIG. 1: Schematic illustration of the jet background subtraction method. 



examine in detail the HIJING fragmentation jet energy 
distribution for those associated with different selected 
fully reconstructed jet energies. In principle, one could 
introduce no such arbitrary definition and put everything 
into a response matrix down to the lowest energy scales. 
In practice, if there are substantial contributions of very 
low energy HIJING fragmentation jet energies to high 
energy reconstruction jets it will be nearly impossible to 
control the systematics and unfold such a matrix. 

Results presented here are based on an analysis of 
750M minimum bias HIJING events. 



IV. RESULTS 

In order to illustrate the background subtraction pro- 
cedure, we show a selection of event displays. Figure [2] 
shows a true dijet pair with R = 0.2 where both jets have 
been matched to reconstructed jets. The reconstructed 
jet has an axis within AR < 0.1 of the true. Also shown 
in the event are the next highest Et reconstructed jet. 
This jet is not matched to any true jets with Et > 5 GeV 
and has Et in the region where we expect fake jets to 
dominate. Figure [3] shows a fake jet with Et =30 GeV 
which is not matched to any true jet from the HIJING 
event. One other fake jet, also not matched to any true 
jets, is shown on the plot. 

We concentrate on central collisions where the underly- 



ing event background is largest. For this study we define 
collision centrality in the HIJING events by the num- 
ber of charged particles with psuedorapidity 3< r\ <4. 
Figure 0] shows the efficiency of finding matches to true 
jets in the most central 10% of collisions for the various 
anti-fcx R parameters as a function of the true jet Et- 
For all R parameters the efficiency rises with jet Et and 
approaches 100% between 20 and 30 GeV. 

However, in order to quantify the jet performance we 
need to understand the contribution to the reconstructed 
jet Et spectrum from jets which are not matched to any 
true HIJING jet, "fake jets" . In Figure [5] we show the 
true, reconstructed and fake jet Et spectra for R = 0.2 
(left), 0.3 (middle) and 0.4 (right) for the 10% most cen- 
tral Au+Au at ^/JWn= 200 GeV HIJING events. Shown 
as red are the true HIJING fragmentation jet distribu- 
tions. The points show the final reconstructed jet dis- 
tribution. This is broken down into those jets that are 
matched with a true HIJING jet and those that are not 
matched with a true HIJING jet. To be considered 
matched the jet axis of the reconstructed jet must be 
within AR < 0.25 of the true HIJING jet and the HI- 
JING jet must have Et >5 GeV. One observes a good 
match between true HIJING and matched reconstructed 
jet distributions taking into account the additional en- 
ergy resolution blurring from the underlying event sub- 
traction. One observes a very large contribution fraction 
of reconstructed jets are not matched at low Et ', the 
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FIG. 2: Event display for a Et = 18 GeV true dijet pair 
jet matched to a Et=15 and 16 GeV reconstructed jets in a 
b=1.8fm HIJING event. All jets shown in this event display 
are reconstructed with the anti-fcr algorithm with R=0.2. 
Both histograms show the background subtracted 0.1x0.1 
r] — <f> tower energy. A minimum Et cut of 5 GeV is placed 
on all jets shown in this display. The stars in the left panel 
show the true HIJING jets and box below shows jet Et and 
rj,(f> location of the jet axis. The right panel shows the recon- 
structed jets. The jet labeled Rl is reconstructed at Et =16 
GeV and matched to the HI jet in the left panel and R2 is 
matched to H2. The other reconstructed jet in the event with 
Et > 5 GeV is shown as R3. It is not associated with any 
true jets with Et > 5 GeV. 



fraction then falls quickly and goes below the matched 
reconstructed jets at around 18 GeV in the R = 0.2 case. 
The crossing point is at higher Et for R=0.3 and 0.4 jets, 
around 25 and 30 GeV, respectively. 

The low fake rate alone is not sufficient due to the ar- 
bitrariness of the 5 GeV separation between "fake" and 
"true" jet fragmentation associations. Shown in Figure[S] 
is the distribution of HIJING true energies for fully re- 
constructed jet energies (with different selections). The 
upper left panel for reconstructed jets with R = 0.2 and 
energies 15-20 GeV shows a peaked distribution around 



FIG. 3: Event display for a E T = 30 GeV fake jet. All jets 
shown in this event display are reconstructed with the anti-fcr 
algorithm with R=0.4. Both histograms show the background 
subtracted 0.1 x 0.1 T] — (f> tower energy. A minimum Et cut 
of 5 GeV is placed on all jets shown in this display. There 
are no true HIJING jets with Et >5 GeV in this event. The 
right panel shows the reconstructed jets. The jet labeled Rl 
is reconstructed at Et = 30 GeV. The other reconstructed 
jet in the event with Et > 5 GeV is shown as R2. The jet 
Et and rj,<f> locations are shown in the bottom right box. 



~ 15 GeV. The tail to lower energies could in princi- 
ple be accounted for in a response matrix (though with 
great care and systematic cross checks). However, as 
one moves to higher energies 25-30 GeV, there is essen- 
tially no tail contributions and a peak around 26 GeV and 
width of 5 GeV. This indicates a regime where a stan- 
dard response matrix and unfolding procedure should be 
successful. Similar plots are shown for R = 0.3 and R = 
0.4. There is a shift downward from the reconstructed 
jet energies to the corresponding true jet energies due to 
the rapid fall off of the jet cross section with energy and 
a tail to low HIJING jet energies that disappears with 
increasing reconstructed jet Et and the corresponding 
decrease of fake jets. 

In order to quantify the purity of the reconstructed jet 
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FIG. 4: Matching efficiency for true HIJING jets as a func- 
tion of jet Et for anti-fcr R parameters 0.2 (black), 0.3 (red) 
and 0.4 (blue). For a jet to be considered matched the recon- 
structed jet axis must be within AT? <0.25 of the HIJING jet 



sample, we have fit the distributions with a background 
contribution which falls exponentially with increasing jet 
Et and a Gaussian with a free mean and width. The 
results of those fits, along with the fractions of the total 
reconstructed jets (both matched and unmatched) which 
are included in the Gaussian are shown in Tables Hl lIIII 



< E T ,reco > (GeV) 


< E T ,true > (GeV) 


OE T 


s 

S+B 


17.2 ± 0.00602 


14.6 ± 0.0537 


3.77 ±0.0306 


0.211 


22.1 ±0.012 


18.6 ±0.051 


4.22 ±0.038 


0.745 


27.2 ±0.0271 


24.1 ±0.0815 


3.79 ± 0.0526 


0.884 


32.3 ± 0.0535 


29.0 ±0.152 


4.21 ±0.121 


0.925 



TABLE I: Jet parameters from fits to the plots in Figure[6]for 
R = 0.2 anti-fcr jets with centrality from - 10%. < ET,reco > 
is the mean reconstructed jet Et within the 5 GeV wide bins. 
< ET,true > and <7ET,true are the mean and width from the 
Gaussian component in the fit and S ^_ B is the fraction of the 
area of the fit that is included in the Gaussian component 
rather than the exponential. 



< E T ,reco > (GeV) 
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22.1 ±0.00748 


14.9 ±0.124 


5.59 ± 0.0634 


0.174 


27.2 ±0.0169 


19.7 ±0.179 


5.83 ±0.144 


0.576 


32.2 ± 0.0354 


23.7 ±0.261 


6.33 ±0.211 


0.802 


37.5 ±0.072 


29.6 ±0.329 


7.04 ± 0.306 


0.908 



TABLE II: Jet parameters from fits to the plots in Figure[7]for 
R = 0.3 anti-fcr jets with centrality from - 10%. < ET,reco > 
is the mean reconstructed jet Et within the 5 GeV wide bins. 
< ET,true > and <7ET,true are the mean and width from the 
Gaussian component in the fit and $ i_ B is the fraction of the 
area of the fit that is included in the Gaussian component 
rather than the exponential. 



In addition to the fake jet contribution to the recon- 
structed jet sample it is also important to quantify the jet 



< E T ,reco > (GeV) 
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<JE T 
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27.1 ±0.01 


15.5 ±0.188 


6.82 ±0.119 


0.128 


32.2 ± 0.0222 


20.3 ±0.412 


7.06 ±0.258 


0.417 


37.3 ± 0.0469 


25.6 ± 0.357 


5.2 ±0.427 


0.379 


42.4 ± 0.0987 


30.4 ±0.616 


4.69 ±0.726 


0.444 



TABLE III: Jet parameters from fits to the plots in Fig- 
ure [H] for R = 0.4 anti-fcr jets with centrality from - 10%. 
< Et.tbco > is the mean reconstructed jet Et within the 5 
GeV wide bins. < ET,true > and <JET,true are the mean and 
width from the Gaussian component in the fit and „^ B is the 
fraction of the area of the fit that is included in the Gaussian 
component rather than the exponential. 



energy resolution and scale for our algorithm. In order to 
do this we have embedded PYTHIA JH] (version 6.421) 
jets into our HIJING events. One PYTHIA event with 
a high pt dijet was embedded into every HIJING event. 
The PYTHIA and reconstructed jets are required to obey 
the same matching cut of AR < 0.25 as the fake jet study 
discussed above. The jet energy resolution and jet energy 
scale are shown in Figure [9] the anti-fey R parameters 0.2 
and 0.4 for central HIJING events and PYTHIA events 
(not embedded into HIJING) put into towers in the same 
manner as the HIJING events. The jet energy resolu- 
tion improves with increasing jet energy and decreasing 
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jet R as expected. The energy scale, 
within ss5% of zero for the anti-fc/r R parameters consid- 
ered here. The energy offset for the PYTHIA jets is due 
to the imposed tower segmentation. For the purposes 
of this study we did not pursue further refinements. A 
similar resolution evaluation has been done by the CMS 
collaboration IJI. 



V. CONCLUSIONS 

We have performed a HIJING study of jet reconstruc- 
tion using an iterative background subtraction method 
and full calorimetric information. We have shown that 
in this case we are able to reconstruct the input HIJING 
jets with a large signal to background with Et > 20 GeV 
for R = 0.2 jets, 30 GeV for R = 0.3 jets and 40 GeV 
for R=0.4 jets. The results presented here are obtained 
without any additionally rejection of fake jets, though it 
is possible the reconstructed jet purities shown here could 
be further increased with fake jet rejection of some kind. 

This study was designed to evaluate the feasibility of 
purely calorimetric jet measurements at RHIC. The re- 
sults here are obtained using an ideal model of the de- 
tector and suggest promise for such measurements. This 
study is of course limited in scope. We have not taken 
into effect any detector effects aside from geometrical ac- 
ceptance and granularity. More detailed studies need to 
be done to demonstrate the suitability of any particular 
detector design for jet measurements. 
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FIG. 5: Et spectra for true HIJING jets (red line) and reconstructed jets (black points). The reconstructed jets are further 
divided into those which are matched to a true HIJING jet (blue line) and those which are not matched to a true HIJING jet 
("fake jets", black line). To be considered matched the axis of the true HIJING jet and the reconstructed jet must be within 
AR <0.25 and the HIJING jet must have Et >5 GeV. Shown are results for 0-10% central HIJING events using anti-fcr jets 
with R=0.2 (a), R=0.3 (b) and R=0.4 (c). 
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FIG. 6: True E T for reconstructed jets anti-fc T R = 0.2 for the reconstructed jet E T 15-20 GeV (a), 20-25 GeV (b), 25-30 GeV 
(c) and 35-40 GeV (d). The lines show the results of fits containing a background component which is exponentially falling 
(dashed line) and a signal Gaussian component (dot-dashed line). The total fit is shown as a solid line. The plots show the „7_ B 
where the signal (S) is determined from the area under the Gaussian within ±2<r of the mean and the total background (B) 
includes both those jets reconstructed with a > 2<r energy mismatched and those which were not matched at all to a HIJING 
jet. Fit parameters are shown in Tabled] 
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FIG. 7: True E T for reconstructed jets anti-fc T R = 0.3 for the reconstructed jet E T 20-25 GeV (a), 25-30 GeV (b), 30-35 GeV 
(c), 35-40 GeV (d). The lines show the results of fits containing a background component which is exponentially falling (dashed 
line) and a signal Gaussian component (dot-dashed line). The total fit is shown as a solid line. The plots show the S ^ B where 
the signal (S) is determined from the area under the Gaussian within ±2a of the mean and the total background (B) includes 
both those jets reconstructed with a > 2<r energy mismatched and those which were not matched at all to a HIJING jet. Fit 
parameters are shown in Table [Til 
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FIG. 9: Jet energy resolution as a function of the true jet energy for jets reconstructed in 0-10% central HIJING events. Shown 
in the plot are anti-fcr R = 0.2 (black), R = 0.3 (red) and R = 0.4 (blue). 



